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ABSTRACT 

The nature of the progenitors of Type la supernovae (SNe la) is still unclear. In this 
paper, by considering the effect of the instability of accretion disk on the evolution 
of white dwarf (WD) binaries, we performed binary evolution calculations for about 
2400 close WD binaries, in which a carbon-oxygen WD accretes material from a main- 
sequence star or a slightly evolved subgiant star (WD + MS channel), or a red-giant 
star (WD -I- RG channel) to increase its mass to the Chandrasekhar (Ch) mass limit. 
According to these calculations, we mapped out the initial parameters for SNe la in 
the orbital period-secondary mass (logP' — Mj) plane for various WD masses for 
these two channels, respectively. We confirm that WDs in the WD -I- MS channel with 
a mass as low as 0.61 Af© can accrete efficiently and reach the Ch limit, while the 
lowest WD mass for the WD -t- RG channel is 1.0 Mq. We have implemented these 
results in a binary population synthesis study to obtain the SN la birthrates and the 
evolution of SN la birthrates with time for both a constant star formation rate and 
a single starburst. We find that the Galactic SN la birthrate from the WD -I- MS 
channel is '^1.8 x 10~^ yr"""^ according to our standard model, which is higher than 
previous results. However, similar to previous studies, the birthrate from the WD + 
RG channel is still low (~3 x 10~^ yr"^)- We also find that about one third of SNe la 
from the WD -I- MS channel and all SNe la from the WD + RG channel can contribute 
to the old populations (>1 Gyr) of SN la progenitors. 
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1 INTRODUCTION 

Type la supernovae (SNe la) play an important role in as- 
trophysics, especially in cosmology. They appear to be good 
cosmological distance indicators and have been applied suc- 
cessfully to the task of determining cosmological parameters 
(e.g. Q and A: Riess et al. 1998; Perlmutter et al. 1999). 
They are also a key part of our understanding of galactic 
chemical evolution owing to the main contribution of iron 
to their host galaxies (e.g. Greggio & Renzini 1983; Mat- 
teucci & Greggio 1986). Despite their importance, several 
key issues related to the nature of their progenitors and the 
physics of the explosion mechanisms are still not well under- 
stood (Hillebrandt & Niemeyer 2000; Ropke & Hillebrandt 
2005; Podsiadlowski et al. 2008; Wang et al. 2008a), and 
no SN la progenitor system has been conclusively identified 
pre-explosion. 

There is a broad agreement that SNe la are thermonu- 
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clear explosions of carbon-oxygen white dwarfs (CO WDs) 
in binaries (see the review of Nomoto, Iwamoto & Kishi- 
moto 1997). Over the last few decades, two competing pro- 
genitor models of SNe la were discussed frequently, i.e. the 
double-degenerate (DD) and single-degenerate (SB) mod- 
els. Of these two models, the SB model (Whelan & Iben 
1973; Nomoto, Thielemann & Yokoi 1984; Fedorova, Tu- 
tukov & Yungelson 2004; Han 2008) is widely accepted at 
present. It is suggested that the BB model, which involves 
the merger of two CO WBs (Iben & Tutukov 1984; Web- 
bink 1984; Han 1998), likely leads to an accretion- induced 
collapse rather than a SN la (Nomoto & Iben 1985; Saio 
& Nomoto 1985; Timmes, Woosley & Taam 1994). For the 
SB model, the companion is probably a main-sequence (MS) 
star or a slightly evolved subgiant star (WB -I- MS channel), 
or a red-giant star (WB + RG channel), or even a He star 
(WB -|- He star channel) (Hachisu, Kato & Nomoto 1996; Li 
& van den Heuvel 1997; Hachisu et al. 1999a; Langer et al. 
2000; Han & Podsiadlowski 2004, 2006; Wang et al. 2009a,b). 
Although the SB model is currently a favorable progenitor 
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model of SNe la, any single channel of the model cannot ac- 
count for the birthrate inferred observationally. Note that, 
some recent observations have indirectly suggested that at 
least some SNe la can be produced by a variety of different 
progenitor systems (e.g. Hansen 2003; Ruiz-Lapuente et al. 
2004; Patat et al. 2007; Voss & Nelemans 2008; Wang et al. 
2008b; Justham et al. 2009). 

At present, various progenitor models of SNe la can be 
examined by comparing the distribution of the delay time 
(between the star formation and SN la explosion) expected 
from a progenitor channel with that of observations (e.g. 
Chen & Li 2007; Meng, Chen & Han 2009; Lii et al. 2009; 
Ruiter, Belczynski & Fryer 2009). Mannucci, Delia VaUe & 
Panagia (2006) argued for the existence of two separate SN 
la populations, a 'prompt' component with a delay time less 
than ~100 Myr, and a 'delayed' component with a delay time 
~3 Gyr. By investigating the star formation history (SFH) 
of 257 SN la host galaxies, Aubourg et al. (2008) found 
evidence of a short-lived population of SN la progenitors 
with lifetimes less than 180 Myr. Botticella et al. (2008) and 
Totani et al. (2008) analyzed host galaxies of SNe la and 
concluded that a substantial fraction of SNe la must have 
long delay times on the order of 2—3 Gyr. Moreover, Schaw- 
inski (2009) recently constrained the minimum delay time of 
21 nearby SNe la by investigating their host galaxies (early- 
type galaxies). The study showed that these early- type host 
galaxies lack 'prompt' SNe la with a delay time less than 
~100 Myr and that ~70 per cent SNe la have minimum de- 
lay times of 275 Myr— 1.25 Gyr, while at least 20 per cent 
SNe la have minimum delay times of at least 1 Gyr at 95 
per cent confidence and two of these four SNe la are likely 
older than 2 Gyr. 

For SNe la with the short delay times, Wang et al. 
(2009a) studied a WD + He star channel to produce SNe 
la, in which a CO WD accretes material from a He MS star 
or a He subgiant to increase its mass to the Chandrasekhar 
(Ch) mass. The study showed the parameter spaces for the 
progenitors of SNe la. By using a detailed binary population 
synthesis (BPS) approach, Wang et al. (2009b, WCMH09) 
found that the Galactic SN la birthrate from this channel is 
~0.3 X 10~^ yr~^, and that this channel can produce SNe 
la with short delay times (~45— 140 Myr). For SNe la with 
the long delay times (>lGyr), this requires that the mass 
of the companion should be <2Mq. 

Recently, Xu & Li (2009) emphasized that the mass- 
transfer through the Roche lobe overflow (RLOF) in the evo- 
lution of WD binaries may become unstable (at least during 
part of the mass-transfer lifetime) , i.e. the mass-transfer rate 
is not equivalent to the mass-accretion rate onto the WD. 
This important feature has been ignored in nearly all of the 
previous theoretical works on SN la progenitors except for 
King, Rolfe & Schenker (2003) Q and Xu & Li (2009), who 
inferred that the mass-accretion rate onto the WD during 
dwarf nova outbursts can be sufficiently high to allow steady 
nuclear burning of the accreted matter and growth of the 
WD mass. In particular, the study of Xu & Li (2009) en- 
larges the region of SN la parameter spaces. However, they 
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only give the SN la parameter spaces with WD initial masses 
of 0.8 and 1 Mq. More detailed work is obviously needed to 
investigate the influence of the accretion disk on the flnal 
results and to give SN la birthrates and delay times. 

Including the effect of the instability of accretion disk on 
the evolution of WD binaries, the purpose of this paper is to 
study the WD -I- MS and WD -I- RG channels towards SNe la 
comprehensively and systematically, and then to determine 
the parameter spaces for SNe la, which can be used in BPS 
studies. In Section 2, we describe the numerical code for 
the binary evolution calculations and the grid of the binary 
models. The binary evolutionary results are shown in Section 
3. We describe the BPS method in Section 4 and present 
the BPS results in section 5. Finally, a discussion is given in 
Section 6, and a summary in Section 7. 



2 BINARY EVOLUTION CALCULATIONS 

In our WD binary models, the lobe-fllling star is a MS star 
or a subgiant star (WD -I- MS channel), or a red-giant star 
(WD -|- RG channel) . The star transfers some of its material 
onto the surface of the WD, which increases the mass of the 
WD as a consequence. If the WD grows to 1.378 Mq, we 
assume that it explodes as a SN la. 



2.1 Stellar evolution code 

We use Eggleton's stellar evolution code (Eggleton 1971, 
1972, 1973) to calculate the WD binary evolutions. The code 
has been updated with the latest input physics over the 
last three decades (Han, Podsiadlowski & Eggleton 1994; 
Pols et al. 1995, 1998). RLOF is treated within the code 
described by Han, Tout & Eggleton (2000). We set the ratio 
of mixing length to local pressure scale height, a = I /Hp, 
to be 2.0 and set the convective overshooting parameter, 
Sov, to be 0.12 (Pols et al. 1997; Schroder, Pols & Eggleton 
1997), which roughly corresponds to an overshooting length 
of ~0.25 pressure scaleheights (Hp). The opacity tables are 
compiled by Chen & Tout (2007) from Iglesias & Rogers 
(1996) and Alexander & Ferguson (1994). In our calculations 
we use a typical Pop I composition with H abundance X = 
0.70, He abundance Y = 0.28 and metaUicity Z = 0.02. 



2.2 Mass-accretion rates 

During the mass-transfer through the RLOF, the accreting 
material can form an accretion disk surrounding the WD, 
which may become thermally unstable when the effective 
temperature in the disk falls below the H ionization temper- 
ature ~6500K (e.g. van Paradijs 1996; King et al. 1997; La- 
sota 2001). This also corresponds to a critical mass-transfer 
rate below which the disk is unstable. Similar to the work of 
Xu & Li (2009), we also set the critical mass-transfer rate 
for a stable accretion disk to be 

Mcr.disk ^ 4.3 X 10"^ (Porb/4hr)i-^M0yr-\ (1) 

for WD accretors (van Paradijs 1996), where Porb is the or- 
bital period. The locations of various types of cataclysmic 
variable stars (e.g. the UX UMa, U Gem, SU UMa, and 
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Z Cam systems) in a (Porb, Afcr.disk) diagram are well de- 
scribed by this expression (Smak 1983; Osaki 1996; van 
Paradijs 1996). 

If the mass-transfer rate, IM2I, is higher than the criti- 
cal value Mcr,disk, wo assume that the accretion disk is stable 
and the WD accretes smoothly at a rate Mace = \M2\; other- 
wise the WD accretes only during outbursts and the mass- 
accretion rate is Mace = |M2|/d, where d is the duty cy- 
cle. The mass-accretion rate is Mace = during quiescence. 
King, Rolfc & Sclicnkcr (2003) showed that for typical val- 
ues of the duty cycle ~0.1 to a few 10~^ the accretion rates 
onto the WD during dwarf nova outbursts can be sufficiently 
high to allow steady nuclear burning of the accreted matter. 
The limits on the duty cycles of dwarf nova outbursts come 
from observations (Waxner 1995): (1) The outburst intervals 
for each object are quasi-periodic, but within the dwaxf nova 
family, the intervals can range from days to decades. (2) The 
lifetime of an outburst is typically from 2 to 20 days and is 
correlated with the outburst interval. The quasi-periodicity 
of the dwarf nova outbursts allows to use a single duty cy- 
cle to roughly describe the change in the mass-transfer rate, 
though we note that this is a simplification of the real, com- 
plicated processes. In this paper we sot the duty cycle to bo 
0.01 (we will discuss the effects of varying the assumed value 
in the discussion section). 

2.3 Mciss-growth rates 

Instead of solving stellar structure equations of a WD, we 
adopt the prescription of Hachisu et al. (1999a) for the mass- 
growth of a CO WD by accretion of H-rich material from 
its companion. The prescription is given below. If the mass- 
accretion rate of the WD, Mace, is above a critical rate, 
Mcr.wD, we assume that the accreted H steadily burns on 
the surface of the WD and that the H-rich material is con- 
verted into He at a rate Mcr.wo- The unprocessed matter 
is assumed to be lost from the system as an optically thick 
wind at a mass- loss rate Mvind = IM2I — Mcr,wD. The crit- 
ical mass-accretion rate is 

Mer.wD = 5.3x10"'^ ^^-'^^ (Mwd/Mq-0.4) Mq yr-\ (2) 

where X is the H mass fraction and Mwd is the mass of the 
accreting WD. 

The following assumptions are adopted when |Mace| is 
smaller than Afcr.wD. 

(1) When |Maec| is less than Mer,wD but higher than 
|Mcr,wD, the H-shell burning is steady and no mass is lost 
from the system. 

(2) When |Afacc| is lower than |Mcr,WD but higher than 
|Mcr,wD, a very weak H-shell flash is triggered but no mass 
is lost from the system. 

(3) When |Maee| is lower than ^M^ ,WD, the H-shell 
flash is so strong that no material is accumulated onto the 
surface of the WD. 

We define the mass-growth rate of the He layer under 
the H-shell burning as 

Mhc = J7H|Macc|, (3) 

where rin is the mass-accumulation efficiency for H-shell 
burning. According to the assumptions above, the values 
of rju are: 



{Mer,WD/|Maee|, |Maec| > Mcr.WD, 
1, Mer.WD > |Maee| > |Mer,WD, (4) 

0, |Maee| < |Mer ,WD. 

When the mass of the He layer reaches a certain value. 
He is assumed to be ignited. If He-shell flashes occur, a part 
of the envelope mass is assumed to be blown off. The mass- 
growth rate of WDs in this case is linearly interpolated from 
a grid computed by Kate & Hachisu (2004), where a wide 
range of WD mass and mass-accretion rate were calculated 
in the He-shell flashes. 

We define the mass-growth rate of the CO WD, Mco, 

as 

Mco = 77HcMhc = »7Ho77H |Mace | , (5) 

where rjue is the mass-accumulation efficiency for He-shell 
flashes. 



2.4 Orbital angular momentum losses 

The evolution of these WD binaries is driven by the nuclear 
evolution of the donor stars, and the change of the orbital 
angular momentum of the binaries is mainly caused by wind 
mass-loss from the WD. We assume that the mass lost from 
these binaries carries away the same specific orbital angular 
momentum of the WD (the mass-loss in the donor's wind 
is supposed to be negligible, but its effect on the change 
of the orbital angular momentum, i.e. magnetic braking, is 
included, e.g. Li & van don Hcuvol 1997). 

For the magnetic braking (MB) effect, we adopt the de- 
scription of angular momentum loss from Sills, Pinsonneault 
& Terndrup (2000), 

— S / X (f.S / s -0.5 (,DJ 

where K = 2.7 x lO*'^ gcm^s (Andronov, Pinsonneault & 
Sills 2003), uj is the angular velocity of the binary, and 
Wcrit is the critical angular velocity at which the angular 
momentum loss rate reaches a saturated state (Krishna- 
murthi et al. 1997). Following the suggestion of Podsiad- 
lowski, Rappaport & Pfahl (2002), we also add an ad hoc 
factor cxp(— 0.02/gconv + 1) if qconv < 0.02 in Equation (6), 
where geonv is the mass fraction of the surface convective en- 
velop, to reduce the MB effect when the convective envelope 
becomes too small. 



2.5 Grid calculations 

We incorporate the prescriptions above into Eggleton's stel- 
lar evolution code and follow the evolution of WD -|- MS 
systems. We have calculated about 2400 WD + MS sys- 
tems, and obtained a large, dense model grid, in which the 
lobe-filling star is a MS star or a subgiant star (WD -|- MS 
channel), or a KG star (WD -I- RG channel). The initial 
mass of the donor star, M2, ranges from 0.5— 4.0 Mq; the 
initial mass of the CO WD, M^^d, is from 0.61-1.20 M©; 
the initial orbital period of the binary system, P\ changes 
from the minimum value, at which a zero-age MS (ZAMS) 
star would fill its Roche lobe, to ~40 d. 
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Figure 1. An example of binary evolution calculations. In panel (a), the solid and dash-dotted curves show M2 and Mwd varing with 
time, respectively. The open circles represent the evolution of Mco during outbursts. In panel (b), the evolutionary track of the donor 
star is shown as a solid curve and the evolution of orbital period is shown as a dash-dotted curve. Dotted vertical lines in both panels 
and asterisks in panel (b) indicate the position where the WD is expected to explode as a SN la. The initial binary parameters and the 
parameters at the moment of the SN la explosion are also given in these two panels. 



3 BINARY EVOLUTION RESULTS 

3.1 An example of binary evolution calculations 

In Fig. 1, we present an example of binary evolution cal- 
culations. Panel (a) shows the M2, Mco and Mwd varing 
with time, while panel (b) is the evolutionary track of the 
donor star in the Hcrtzsprung-Russcll diagram, where the 
evolution of the orbital period is also shown. The binary 
shown in this case is (M2, MJ^q, log(PYday)) = (1.6, 1.0, 
—0.1), where M2, M^d and P' are the initial mass of the 
donor star and the CO WD in solar masses, and the initial 
orbital period in days, respectively. The donor star fills its 
Roche lobe on the MS which results in case A mass-transfer. 
During the whole evolution, the mass-transfer rate is lower 
than the critical value given by Equation (1). Thus, the ac- 
cretion disk experiences instability. The mass-accretion rate 
Mace of the WD exceeds ^Mcr.wD after the onset of RLOF, 
leading to the mass-growth of the WD. With the contirmous 
decreasing of \M2\, when Mace drops below |Mcr,wD after 
about 5 X 10^ yr, the H-shell flash is so strong that no ma- 
terial is accumulated onto the surface of the WD. At about 
1.7 X 10^ yr, [Mai increases again as the donor star evolves 
off the MS and expands, allowing the mass-growth of the 
WD during outbursts. When the WD grows to 1.378 M©, 
it explodes as a SN la. At the SN explosion moment, the 
mass of the donor star is Mf* = 0.6253 M© and the orbital 
period log(P^^/day) = 0.2289. 

3.2 Initial parameters for SN la progenitors 

Figs 2-5 show the final outcomes of about 2400 binary evolu- 
tion calculations in the initial orbital period-secondary mass 
(log P' — M2) plane. The filled squares, circles and stars de- 
note that the WD explodes as a SN la in the optically thick 
wind phase, in the stable H-shell burning phase and in the 
weak H-shell flash phase, respectively. The filled triangles 
represent that the WD explodes as a SN la, in which the 
accretion disk experiences instability. Some binaries fail to 
produce SNe la owing to nova explosions (which prevents 
the WD growing in mass; the open circles in these figures) 



or dynamically unstable mass-transfer (resulting in a com- 
mon envelope; the crosses in these figures). 

The contours of initial parameters for producing SNe la 
arc also presented in these figures (solid curves are for the 
contours from the WD -|- MS channel, while dotted curves 
from the WD -|- RG channel). The loft boundaries of the 
contours in these figures (solid lines; Figs 2-4 and panel (d) 
of Fig. 5) arc set by the condition that RLOF starts when the 
donor star is on the ZAMS, while systems beyond the right 
boundary experience mass-transfer at a very high rate due 
to the rapid expansion of the donor stars in the Hertzsprung 
gap (HG) (solid lines) or RG stage (dotted lines) and they 
lose too much mass via the optically thick wind, preventing 
the WDs increasing their masses to the Ch mass. The upper 
boundaries are set mainly by a high mass-transfer rate owing 
to a large mass-ratio. The lower boundaries arc constrained 
by the facts that the mass-transfer rate Ad-z should be high 
enough to ensure that the WD can grow in mass during 
outbursts and that the donor should be sufficiently massive 
for enough mass to be transferred onto the WD. 

In Fig. 6, we overlay the contours for SN la production 
in the (logP', M2) plane for difi^erent initial WD masses. In 
panel (a), wo give the contours of the WD -I- MS channel 
for producing SNe la with various WD masses (i.e. M^n = 
0.61, 0.65, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2 M©). Note that the 
enclosed region almost vanishes for Af^^ = 0.61 Mq, which 
is then assumed to be the minimum WD mass for producing 
SNe la from this channel. As a comparison, we also give the 
contours of the WD -|- RG channel in panel (b). If the initial 
parameters of a WD binary are located in the contours, a 
SN la is then assumed to be produced. Thus, these contours 
can be expediently used in BPS studies. The data points and 
the interpolation FORTRAN code for these contours can be 
supplied on request by contacting BW. 

4 BINARY POPULATION SYNTHESIS 

In order to investigate SN la birthrates and delay times for 
the WD + MS channel, we have performed a series of Monte 



Progenitors of SNe la with long delay times 5 





X X 


X 


1 






X 




1 






X 




1 

X 




X 


1 . . . . 1 . . . 




X X 


X 






X 


























X X 






X 


X 










■-A 


A 


i\ 






X 


■ ; SN Igl, opticBilly Uiiclc wind 




















if 




A 


A ^ 










X X 




















A 


A 








• ! SN Is, st&bls H — shell burning 




X X 


X 




















A 




































A 










•k' f^N Ta wpfllr H— shpll flash 




X X 




























X 


X 






X 




























X 


X 


SN la, disk instability 
































X 


X 






























/x 






; no SN I& nova — 




























X 


X 


X 






























X 


X 


X 


X ; no SN I&, dyn&micEtlly unst&bls 


























' X 


X 


X 


X 




























X 


X 


X 


X 






























X 






























X 


X 


X 


X 


X 






























X 
























A 


io 


X 


X 


X 


X 


X 


X 






★ 


★ 


• 


• 


★ 


A 


A 


a/ 


X 


X 


X 


X 


X 


X 


X 






★ 


★ 


★ 


★ 


A 


A 


A 




X 


X 


X 


X 


X 


X 


X 






A 


A 


★ 


A 


A 


A 


A , 








X 




X 


X 


X 






A 


A 


A 


A 


A 


A 
















X 


X 






A 


A 


A 


A 


A 


A , 


















X 


X 


































X 


X 






\ A 


A 


A 


A 




a' 




-x. 


X 












X 








A 


A 


Ay 


^A 


A 


A 


A 












A_^ 










A 


A 




A 


A 


A 


A 


A 


A 


A 


A 


A 


A 


*♦ 












A 


A 


A 


A 


A 


A 


A 


A 


A 


A 


A 


A* 








































- 













































b • 










































A 


A 


A 


A 


A 


A 


A 


A 


A 


















*A 


. A 


A 


A 


A 


A 


A 


A 


A 
























•Am 


• A 


A • 


■A. 


•Ap 


• A - 


i. a ° o 








1 










1 










1 


O 


o 


o o o o 





0.5 1 1.5 2 



log (days) 

Figure 2. Final outcomes of the binary evolution calculations in the initial orbital period-secondary mass (logP', Mj) plane of the CO 
WD + MS system for an initial WD mass of 1.2 Mq. Solid curves are for the contours from the WD + MS channel, while dotted curves 
from the WD + RG channel. The filled squares, circles and stars denote that the WD explodes as a SN la in the optically thick wind 
phase, in the stable H-shell burning phase and in the weak H-shell flash phase, respectively. The filled triangles represent that the WD 
explodes as a SN la, in which the accretion disk experiences instability. Open circles indicate systems that experience novae, preventing 
the WD from reaching 1.378 Mq, and crosses are those under dynamically unstable mass-transfer. 
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Figure 3. Similar to Fig. 2, but for an initial WD mass of 1.1 Mq. 
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Figure 4. Similar to Fig. 2, but for initial WD masses of 0.9 and 1.0 Mq. 
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Figure 6. Panel (a): regions in the initial orbital period-secondary mass plane (logP', M2) for the WD + MS channel that produce 
SNe la for initial WD masses of 0.61, 0.65, 0.7, 0.8, 0.9, 1.0, 1.1 and 1.2 M©. The region almost vanishes for = 0.61 M©. Panel 

(b): similar to the panel (a), but the regions are for the WD + RG channel. 
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Carlo simulations in the BPS study. In each simulation, by 
using the Hurley's rapid binary evolution code (Hurley, Pols 
& Tout 2000, 2002), we have followed the evolution of 1 x lO'^ 
sample binaries from the star formation to the formation 
of the WD + MS systems according to three evolutionary 
channels (Sect. 4.2). We assumed that, if the parameters of a 
CO WD + MS system at the onset of the RLOF are located 
in the SN la production regions (panel (a) of Fig. 6), a SN 
la is produced. Note that, the method of the BPS study for 
the WD + RG channel is similar to that of the WD + MS 
channel. 



4.1 Common envelope in binary evolution 

In the SD model of SNe la, the progenitor of a SN la is 
a close WD binary system, which has most likely emerged 
from the common envelope (CE) evolution (Paczyriski 1976) 
of a giant binary system. The CE ejection is still an open 
problem. Similar to the work of WCMH09, we also use the 
standard energy equations (Webbink 1984) to calculate the 
output of the CE phase. The CE is ejected if 



2af 



2ai 



(7) 



where A is a structure parameter that depends on the evo- 
lutionary stage of the donor, Mdon is the mass of the donor, 
A^acc is the mass of the accretor, a is the orbital separation, 
Menv is the mass of the donor's envelope, Rdon is the ra- 
dius of the donor, and the indices i and f denote the initial 
and final values, respectively. The right side of the equation 
represents the binding energy of the CE, the left side shows 
the difference between the final and initial orbital energy, 
and ace is the CE ejection efficiency, i.e. the fraction of the 
released orbital energy used to eject the CE. For this pre- 
scription of the CE ejection, there are two highly uncertain 
parameters (i.e. ace and A). As in previous studies, we com- 
bine ace and A into one free parameter aceA, and set it to 
be 0.5 and 1.5 (e.g. WCMH09). 



4.2 Evolutionary channels to WD -|- MS and WD 

-|- RG systems 

According to the evolutionary phase of the primordial pri- 
mary at the beginning of the first RLOF, there are three 
channels which can form CO WD -|- MS systems and then 
produce SNe la. 

(1) He star channel. The primordial primary first fills 
its Roche lobe when it is in the HG (so-called early Case 
B mass-transfer defined by Kippenhahn & Weigert 1967). 
A CE may be formed when the lobe-filling star evolves to 
the RG stage owing to a large mass-ratio or a convective 
envelope of the mass donor star. After the CE ejection, the 
primary becomes a He star and continues to evolve. After 
the exhaustion of central He, the He star which now contains 
a CO-core may fill its Roche lobe again due to expansion of 
the He star itself, and transfer its remaining He-rich enve- 
lope to the MS companion star, eventually leading to the 
formation of a CO WD -|- MS system. For this channel, 
SN la explosions occur for the ranges Mi,i ~ 4.0— 7.0 Mq, 
M2,i ~ 1.0-2.0 Mq, and ~ 5-30days, where Mi,i, M2,i 



and P' are the initial mass of the primary and the secondary 
at ZAMS, and the initial orbital period of a binary system. 

(2) EAGB channel. If the primordial primary is on 
the early asymptotic giant branch (EAGB, i.e. He is ex- 
hausted in the centre of the star while thermal pulses have 
not yet started), a CE will be formed because of dynam- 
ically unstable mass-transfer. After the CE is ejected, the 
orbit decays and the primordial primaxy becomes a He 
RG. The He RG may fill its Roche lobe and start mass- 
transfer, which is likely stable and leaves a CO WD + 
MS system. For this channel, SN la explosions occur for 
the ranges Mi,i ~ 2.5 - 6.5 Mq, M2,i ~ 1.5 - 3.0 M© and 
P' ~ 200 - 900 days. 

(3) TPAGB channel. The primordial primary fills its 
Roche lobe at the thermal pulsing asymptotic giant branch 
(TPAGB) stage. A CE is easily formed owing to dynam- 
ically unstable mass-transfer during the RLOF. After the 
CE ejection, the primordial primary becomes a CO WD, 
then a CO WD + MS system is produced. For this channel, 
SN la explosions occur for the ranges Mi,i ~ 4.5— 6.5 Mq, 
M2,i ~ 1.5-3.5 Mq, and P' > 1000 days. 

Additionally, there is one channel which can form CO 
WD + RG systems and then produce SNe la (i.e. the 
TPAGB channel above). After a CO WD -I- MS system 
is produced via the TPAGB channel, the MS companion 
star continues to evolve until the RG stage, i.e. a CO WD 
-|- RG system is formed. For the CO WD -|- RG systems, 
SN la explosions occur for the ranges Mi,i ~ 5.0— 6.5 Mq, 
M2,i - 1.0-1.5 Mq, and P' > 1500 days. 

4.3 Basic parameters for Monte Carlo simulations 

In the BPS study, the Monte Carlo simulation requires as in- 
put the initial mass function (IMF) of the primary, the mass- 
ratio distribution, the distribution of initial orbital separa- 
tions, the eccentricity distribution of binary orbit, and the 
star formation rate (SFR) (e.g. Han et al. 2002, 2003; Han, 
Podsiadlowski & Lynas-Gray 2007; WCMH09). 

(1) The IMF of Miller & Scalo (1979) is adopted. The 
primordial primary is generated according to the formula of 
Eggleton, Tout & Fitechett (1989), 

MP = rs^ 

1 (l-JC)0-75 + 0.032(1 -A:)0-25' 

where X is a random number uniformly distributed in the 
range [0, 1] and Mf is the mass of the primordial primary, 
which ranges from 0.1 Mq to 100 Mq. The studies of the 
IMF by Kroupa, Tout & Gilmore (1993) and Zoccah et al. 
(2000) support this IMF. 

(2) The initial mass-ratio distribution of the binaries, 
q , is quite uncertain for binary evolution. For simplicity, we 
take a constant mass-ratio distribution (Mazeh et al. 1992; 
Goldberg & Mazeh 1994), 



n{q') = 1, 



< g' < 1, 



(9) 



where q = Mj/Mf. This constant mass-ratio distribution 
is supported by the study of Shatsky & Tokovinin (2002). As 
an alternative mass-ratio distribution we also consider un- 
corrclated binary components, i.e. both binary components 
are chosen randomly and independently from the same IMF 
(uncorrelated). 

(3) We assume that all stars are members of binaries 
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Table 1. Galactic birthrates of SNe la for different simulation 
sets, where sets 1 and 4 are our standard models for the WD + 

MS and WD -I- RO rhanncls. rosT)ccti\'oly. 



Set 


Channel 


Ctce^ 


n{q') 


u (10-3 yr-l) 


1 


WD + MS 


0.5 


Constant 


1.793 


2 


WD + MS 


1.5 


Constant 


1.415 


3 


WD + MS 


0.5 


Uncorrclated 


0.275 


4 


WD + RG 


0.5 


Constant 


0.028 


5 


WD + RG 


1.5 


Constant 


0.026 



and that the distribution of separations is constant in log a 
for wide binaries, where a is separation and falls off smoothly 
at small separation 



• n{a) = 



>(a/oo)" 



a ^ ao, 
ao < a < ai, 



(10) 



and m w 1.2. This distribution implies that the numbers of 

wide binaries per logarithmic interval are equal, and that 
about 50 per cent of stellar systems have orbital periods less 
than 100 yr (Han, Podsiadlowski & Eggleton 1995). 

(4) A circular orbit is assumed for all binaries. The or- 
bits of semidetached binaries are generally circularized by 
the tidal force on a timescale which is much smaller than 
the nuclear timescale. Also, a binary is expected to become 
circularized during the RLOF. 

(5) We simply assume a constant SFR over the last 
15 Gyr or, alternatively, as a delta function, i.e. a single star- 
burst. In the case of the constant SFR, we calibrate the SFR 
by assuming that one binary with a primary more mgissive 
than 0.8 Mq is formed annually (see Iben & Tutukov 1984; 
Han, Podsiadlowski & Eggleton 1995; Hurley, Pols & Tout 
2002). From this calibration, we can get SFR = 5MQyr^"'^ 
(e.g. Willems & Kolb 2004), which successfully reproduces 
the ^®A1 1.809-MeV gamma^ray line and the core-collapse 
SN rate in the Galaxy (Timmes, Diehl & Hartmann 1997). 
It is obvious that for a galaxy with a different SFR the SN 
la rates will be scaled. For the case of the single starburst, 
we assume a burst producing 10^^ Mq in stars, which is a 
typical mass of a galaxy. In fact, a galaxy may have a com- 
plicated star formation history. We only choose these two 
extremes for simplicity. A constant SFR is similar to the 
situation of spiral galaxies (Yungelson & Livio 1998; Han 
& Podsiadlowski 2004), while a delta function to that of 
elliptical galaxies or globular clusters. 



5 THE RESULTS OF BINARY POPULATION 
SYNTHESIS 

We performed five sets of simulations (see Table 1) with 
metallicity Z = 0.02 to systematically investigate Galactic 
birthrates of SNe la, where sets 1 and 4 are our standard 
models for the WD -I- MS and WD + RG channels, respec- 
tively, and with the best choice of model parameters (e.g. 
WCMH09) . We vary the model parameters in the other sets 
to examine their influences on the final results. 

In Fig. 7, we show the evolution of Galactic SN la 
birthrate by adopting Z = 0.02 and SFR = 5 Mgyr-^ both 
for the WD -|- MS channel (thick curves) and the WD -|- 



& 2 



o 

a 

a* 




log(i/yr) 

Figure 7. Evolution of Galactic birthrates of SNe la for a con- 
stant star formation rate (Z = 0.02, SFR = 5M0yr-i). The key 
to the line-styles representing different sets is given in the upper 
left corner. 



RG channel (thin curves). The simulation for the WD -I- MS 
channel gives Galactic SN la birthrate of ~1.8 x 10-^ yr-^ 
according to our standard model (set 1). The result is higher 
by a factor of 2—3 than previous studies by Han & Pod- 
siadlowski (2004, HP04) (~0.6-0.8 x 10-^ yr-^). This is 
mainly attributed to the effect of the accretion disk instabil- 
ity, which can increase the accretion rate onto the WD dur- 
ing outbursts, leading to higher effective WD growth rates 
and a larger fraction of systems becoming SN la progenitors. 
However, the result from this work is still lower than that in- 
ferred observationally (i.e. 3—4 x 10"^ yr~^: van den Bergh 
& Tammann 1991; Cappellaro & Turatto 1997). The simula- 
tion for a constant mass-ratio distribution with ctccA = 1.5 
(set 2) gives SN la birthrate of ~1.4 x 10^"* yi'"\ which 
is lower than the case of accA — 0.5 (the binaries resulted 
from the CE ejections tend to have slightly closer orbits for 
aceA = 0.5 and are more likely to locate in the SN la pro- 
duction region). If we adopt a mass-ratio distribution with 
uncorrelated binary components and OceA = 0.5 (set 3), the 
SN la birthrate from the WD + MS channel will be lower 
by an order of magnitude (see dotted curve in Fig. 7). This 
is because most of the donors in the WD -|- MS channel are 
not very massive which has the consequence that WDs can- 
not accrete enough mass to reach the Ch mass. Similar to 
previous studies (e.g. Yungelson & Livio 1998; HP04), the 
Galactic SN la birthrate from the WD -|- RG channel is still 
low (~3 x IQ-^ yr~^)- The studies in this paper imply that 
the WD + MS and WD + RG channels arc only subclasses 
of SN la production, and there may be some other chan- 
nels or mechanisms also contributing to SNe la, e.g. WD + 
He star channel or double-degenerate channel. Especially, as 
mentioned by WCMH09, the WD + He star channel can give 



a Galactic birthrate of ^^0.3 x 10-^ yr-^, and is considered 
to be an important channel to produce SNe la with short 
delay times. The SN la birthrates shown in this figure seem 
to be so completely flat after the flrst rise, for the specific 
reason see Sect. 3 of WCMH09. 

Fig. 8 displays the evolution of SN la birthrates for a 
single starburst with a total mass of 10^^ Mq. In the fig- 
ure we see that SNe la from the WD -|- MS channel (thick 
curves) have the delay times of ~250 Myr— 6.3 Gyr (while 
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Figure 8. Similar to Fig. 7, but for a single starburst with a total 
mass of 10" M0. 
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Figure 9. Same as Fig. 8, but for different SD models. The SD 
models are from HP04 (dashed curve), WCMH09 (dotted curve) 
and this work (dot-dashed curve is for the WD + MS channel, 
while solid curve is for the WD + RG channel) . All curves are for 
a constant mass-ratio distribution and with OceA = 0.5. 



SNe la from the WD -|- RG channel have the delay times 
of >3Gyr; thin curves). This figure also shows that a high 
value of OceA leads to a systematically later explosion time 
for these two channels (e.g. the cases of OcoA — 0.5 and 
1.5 for a constant mass-ratio distribution). This is because 
a high value of OceA leads to wider WD binaries, and, as 
a consequence, it takes a longer time for the companion to 
evolve to fill its Roche lobe. 

In Fig. 9, we compare the theoretical delay time distri- 
butions (DTDs) of SNe la for different SD models. The SD 
models are from HP04 (dashed curve), WCMH09 (dotted 
curve) and this work (dot-dashed curve is for the WD + 
MS channel, while solid curve is for the WD -I- RG chan- 
nel). In this figure, the WD -I- MS channel of HP04 may 
contribute to the SNe la with intermediate delay times 
(100 Myr— 1 Gyr: for a detailed discussion see Schawinski 
2009). The results of the WD + He star channel from 
WCMH09 may provide a possible ways to explain the SNe 
la with short delay times (< 100 Myr). In addition, for SNe 
la with the long delay times (>lGyr), about one third of 



SNe la from the WD + MS channel (dot-dashed curve) and 
all SNe la from the WD + RG channel (solid curve) can con- 
tribute to the old populations of SN la progenitors. Note, 
Chen & Li (2007) studied the WD + MS channel by consid- 
ering a circumbinary disk which extracts the orbital angu- 
lar momentum from the binary through tidal torques. The 
study cam also provide a possible ways of producing such old 
SNe la 3Gyr). Finally, we emphasize that SNe la may 
be from several different progenitor systems. 



6 DISCUSSION 

The regions (Fig. 6 in this paper) for producing SNe la de- 
pend on many uncertain input parameters, in particular for 
the duty cycle which is poorly known. The main uncertain- 
ties lie in the facts that it varies from one binary system to 
another and may evolve with the orbital periods and mass- 
transfer rates (e.g. Lcisota 2001; Xu & Li 2009). This is the 
reason why we choose an intermediate value (0.01) rather 
than other extreme ones (e.g. 0.1 or 10~^). However, we 
also did some tests for a higher or lower value of the duty 
cycle. We find that the variation of the duty cycle value will 
influence the regions for producing SNe la, especially for the 
WD + RG channel (for a high value (0.1), the regions will be 
shifted to higher period; for a low value (10~^), the regions 
will be shifted to lower period) . For these two extreme duty 
cycle values, the SN la birthrate from the WD -I- RG chan- 
nel will be lower due to the smaller regions for producing 
SNe la. 

In our binary evolution calculations we have not consid- 
ered the influence of rotation on the H-accreting WDs. The 
calculations of Yoon, Langer & Scheithauer (2004) showed 
that the He-shell burning is much less violent when rotation 
is taken into account. This may significantly increase the 
He-accretion efficiency {'que in this paper). Meanwhile, the 
maximum stable mass of a rotating WD may be above the 
standard Ch mass (i.e. the super-Ch mass model: Uenishi, 
Nomoto & Hachisu 2003; Yoon & Langer 2005; Chen & Li 
2009). However, we mainly focus on the standard Ch mass 
explosions of the accreting WDs in this work. 

In our BPS studies we assume that all stars are in bi- 
naries and about 50 per cent of stellar systems have orbital 
periods less than 100 yr. In fact, this is known to be a simpli- 
fication, and the binary fractions may depend on metallicity, 
environment, spectral type, etc. If we adopt 40 per cent of 
stellar systems have orbital periods below 100 yr by adjust- 
ing the parameters in Equation (10), we estimate that the 
Galactic SN la birthrate from the WD + MS channel will 
decrease to be ~1.4 x 10^"^ yr^^ according to our standard 
model (the birthrate from the WD -|- RG channel will de- 
crease to be ~2 X 10~^ yr~^). In addition, Umeda et al. 
(1999) concluded that the upper limit mass of CO cores 
born in binaries is about 1.07 M0. If this value is adopted 
as the upper limit of the CO WD, the SN la birthrate from 
the WD -I- MS channel wiU decrease to be ~1.7 x 10"^ yr"^ 
(the birthrate from the WD -|- RG channel will decrease to 
be ~1 X lO"'^ yr"^). 

The Galactic SN la birthrate from the WD + RG chan- 
nel is ~3 x 10~^ yr~^ according to our standard model, which 
is low compared with observations, i.e. SNe la from this 
channel may be rare. However, further study on this chan- 
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nel is necessary, since this channel may explain sonic SNc 
la with long delay times. In addition, it is suggested that, 
RS Oph and T CrB, both recurrent novae are probable SN 
la progenitors and belong to the WD + RG channel (e.g. 
Belczynski & Mikolajowska 1998; Hachisu, Kato & Nomoto 
1999b; Sokoloski et al. 2006; Hachisu, Kato & Luna 2007). 
Meanwhile, by detecting Na I absorption lines with low ex- 
pansion velocities, Patat et al. (2007) suggested that the 
companion of the progenitor of SN 2006X may be an early 
RG star. Additionally, Voss & Nclomans (2008) studied the 
pre-explosion archival X-ray images at the position of the 
recent SN 2007on, and they considered that its progenitor 
may be a WD + RG system. 



7 SUMMARY 

Employing Egglcton's stellar evolution code with the pre- 
scription of Hachisu et al. (1999a) for the mass-accretion of 
CO WDs, and including the effect of the instability of ac- 
cretion disk on the evolution of WD binaries, we performed 
binary evolution calculations for about 2400 close WD bina- 
ries. The calculated results further confirm that the disk in- 
stability could substantially increase the mass-accumulation 
efficiency for accreting WDs, and cause the possible SNe la 
to occur in systems with <2 Mq donor stars (see also King, 
Rolfe & Schenker 2003; Xu & Li 2009). We find that the 
Galactic SN la birthrate from the WD -|- MS channel is 
~1.8 X 10~^ yr~^ according to our standard model, which 
is higher than previous results. However, similar to previous 
studies, the birthrate from the WD + RG channel is still 
low (^3 X 10^'' yr^^). We also find that about one third 
of SNe la from the WD -I- MS channel and all SNe la from 
the WD -|- RG channel can contribute to the old popula- 
tions (>lGyr) of SN la progenitors. The companion stars 
of SNe la with long delay times in this work would survive 
in the SN explosion and show distinguishing properties. In 
future investigations, we will explore the properties of the 
companion stars after SN explosion, which could be verified 
by future observations. 
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